clear
rng(20)

simulated_data = [];
outer_counter = 1;
repeats = 10;
num_comps = 10;
parpool('local',num_comps)

%%% UNCOMMENT FIRST ROW FOR KAPPA = 0.94 AND SECOND ROW FOR KAPPA = 0.01
params = [0.1207	0.0296	0.1584	1.538	7.108	0.938];
%params = [0.1653 0.0344 0.2467 1.538 6.9405 0.01];
iteration_data = zeros(4,2);
counter=1;
while counter <= repeats
parfor computer_num = 1:num_comps
temp_data = steady_state_simulation_plots(params);
iteration_data = iteration_data+temp_data;
end
counter = counter+1
end
iteration_data = iteration_data./(repeats*num_comps);
simulated_data = [simulated_data iteration_data];


growth_rate_exports_domestic_contracting_exporters = simulated_data(1,1)/simulated_data(2,1)
growth_rate_exports_domestic_expanding_exporters = simulated_data(3,1)/simulated_data(4,1)